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^rr . In the Brandenberger-Vafa scenario of string gas cosmology, the Universe starts as a small torus 

^^ ! of string length dimension filled with a hot gas of strings. In such extreme conditions, in addition 
to the departure from Einstein gravity which is due to the dilaton, one expects higher curvature 
corrections to be relevant. Motivated by this fact, we study the effect of the leading a'^ corrections 
flj ' in type IIB string theory for this scenario. Within the assumptions of: weak coupling, adiabatic 

T^ • evolution and thermodynamical equilibrium, we perturbatively solved the corresponding equations 
of motion in two different cases: (i) the isotropic case which is governed by a single scale factor 

^ . and (ii) the anisotropic case given by two different scale factors. In the first case, we consider two 
regimes (ia) The Hagedorn regime where the string gas equation of state is that of pressureless dust, 
and (ib) the radiation regime. In the second case, (ii), we only considered a radiation-like equation 
of state. We found that the inclusion of a' corrections affects the scale factor(s) in opposite way in 
the Hagedorn and in the radiation regimes, acting as a driving force for the first one and a damping 
force for the second one. This effect is small for reasonable initial conditions and it is only observed 
at early times. Morever it is bigger in the Hagedorn regime than in the radiation regime. We also 
analyzed the fixed points of the system. We found that there exists a stable dS fixed point, which 
does not exist when the corrections are neglected. 
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1 Introduction 

String theory remains the most promising and theoreticahy appeahng candidate for a complete and 
unified description of all Nature's forces including gravity. Given the energy scales where string 
theory is valid, cosmology offers an invaluable, and maybe the only, testing ground for string theory 
ideas. String Gas Cosmology is one of the existing approaches trying to obtain our observed 4D 
universe starting from a string theory set up. This approach was initiated by the pioneering work 
of Brandenberger and Vafa [1]. It aims to address two important questions which are most relevant 
for cosmology: the initial cosmological singularity and the dimensionality of space-time. 

The initial singularity arises from the fact that physical quantities, such as the Ricci scalar and 
temperature, become singular as one extrapolates standard cosmology to very early times. String 
theory has a built-in symmetry called T-duality, which acts as a{t) -^ l/a{t), where a(t) is the 
scale factor, that prevents such a singularity. This means that a Universe with a tiny scale factor 
is equivalent to one with a very large scale factor. In this way, all the physical quantities remain 
well-behaved at all times. This feature is also used in the so-called pre big bang scenario [2] . 

Another important issue is the dimensionality of the space-time. String theory (in its supersym- 
metric version) requires for consistency that the dimension of space-time is 10. If string theory is 
a valid description at early times, then what is the mechanism that made our 3 observed spatial 
dimensions grow and the remaining 6 become unobservably small? In the Brandenberger and Vafa 
scenario all nine spatial dimensions start out toroidally compactified with a common radius of di- 
mension of the string length £s = a'^'^, where a' is the Regge slope related to the string tension 
through T = l/2TTa' . The strings in this background will have both momentum (KK) and winding 
modes (W) whose dynamics will determine how many dimensions expand and how many stay at 
the string length. A crucial feature is that P^-modes impede the expansion of the dimensions they 
are wrapping. The qualitative argument is that it is easier for the W and W to annihilate in the 
maximum dimension of space-time D = 4 because of the naive numerological argument 2-|-2=4. 

Even though the Brandenberger- Vafa (BV) scenario is appealing, there exist obstacles to its 
concrete realization. The dynamics described above presupposes thermal equilibrium, which means 
that the annihilation rate of WW in the large dimensions has to be fast enough compared with the 
expansion of the Universe. Many studies have been done to see if this qualitative behavior survives, 
e.g. in [3,4]. According to [4] there are two possible situations: either all the dimensions remain 
small, or they expand and become large. In the first case, the WW annihilation is not efficient and 
a substantial number of these modes remain preventing the large dimensions to grow. The reason 
behind that is not the fact that the expansion of the Universe is very fast at the time of annihilation, 
but it is instead that the string coupling (which enters the WW annihilation cross section) becomes 
very small in a finite amount of time. In the second case, the universe starts with few M^-modes, 
they annihilate completely and all dimensions decompactify. Only very special initial conditions 
lead to the desired outcome. 

So far, many of the existing studies of this scenario have been done using the low energy effective 
action of type IIB theory - dilaton gravity system neglecting a' corrections [4-36]. However, since 
the dynamics are supposed to take place when the scale factors are of order of the string length, 
it is then important to study the leading a' corrections in the effective action. Due to the length 
scale we are dealing with, it is obvious that all corrections will contribute in the same footing to the 
dynamics. In this paper, and as a modest step towards a full understanding of these corrections, we 



will address the effect of the leading corrections in type II string theory which are proportional to 

This paper is organized as follows: In section [2 we summarize some aspects of the BV scenario 
in the absence of a' corrections. Then we summarize the dynamics for two cases: a) the isotropic 
case, a single scale factor, where we consider two different regimes (equations of state), the Hagedorn 
and the radiation regimes; b) the anisotropic case, with two scale factors, where we consider the 
radiation regime. We analyse the dynamics of the system in the presence of the leading a' corrections, 
providing dilaton gravity equations for the isotropic case in section 3.1. In sections 3.1.1 and 3.1.2, 
we find the perturbative solutions for the Hagedorn and the radiation regime. The fixed points for 
the exact equations of motion and their stability are analysed in sections 3.1.3 and 3.1.4. For the 
anisotropic case, the perturbative solutions are found in section 3.2 and the presence of fixed points 
and their stablity are studied in sections 3.2.1 and 3.2.2. We present our conclusions in the final 
section. In the Appendix, we collect some cumbersome formulae. 

2 The BV scenario without a'-corrections 

In the standard BV scenario the Universe starts out at a very high temperature, close to the Hagedorn 
temperature Th = l/27rv2c7, which is the limiting temperature in type IIB string theory [37]. The 
geometry of the Universe is taken as a 9 dimensional small torus described by the metric 
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ds^ = Gab dx^ dx^ = -dt^ + l^Yl "?(*) "^^l 0<9,<27r (1) 

The universe is assumed to be homogeneous and in thermal equilibrium. The scale factors are 
parametrized as aj(t) = e '^^y'') . The torus radii are assumed to be all more or less equal to a common 
value given by the string length Ig . The Universe is assumed to be filled with a gas of closed strings^ . 
For simplicity, the string gas is assumed to be an ideal gas of strings, whose free energy F{Xi, (3) at 
finite temperature T = (5~^ is computed from the string partition function Z(/3, Aj) (see e.g. [14]) 
using the thermodynamical identity -F(/?, Aj) = — /3^^ logZ(/3, Aj). For our purposes it is enough to 
include only massless string modes. 

To summarize, the standard BV scenario is based on two ingredients: the first ingredient, which 
comes from string theory, is the low energy effective action (at tree level in a' and without form 
fluxes) for type IIB string theory. The second ingredient, which comes from string thermodynamics, 
is the free energy -F(Aj, j3) of the string gas. 

Thus, the action for the BV scenario is 

5b V = Sq + Sm, (2) 

^ / d^Ox ^/^e -^-^ {R + ^{d4>f) , (3) 

fdty^^Q^F{X^,(3y^^G^), (4) 
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^Even in the presence of D-branes of various dimensionalities, the dynamics is still dominated by funda- 
mental strings [7]. 



where (/) = (f>{t) is the dilaton field and kiq is related to the string length through k^q = i(27r)'^a'^. 
The equations governing the dynamics are obtained by varying the action Sby [5] 
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where the shifted dilaton is defined as -0 = 2</> — > Aj. The pressure and the total energy of the 
string gas P and E are defined as 



j=i 
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The conservation of the total energy reads 

9 

E + Y^ XiPi = 0. (10) 

i=l 

Eq. (|1U|) is equivalent to the adiabaticity condition which means that at any given temperature, 
the system will adjust its volume so to keep its entropy constant. Therefore, during the adiabatic 
evolution, there exists a one-to-one relationship {Aj} <-> (5~^ = T. 

It is worth mentioning that Eqs. (|S]|7j) are invariant under T-duality: 
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Aj — > — Aj and (j!) — > — ^^ Aj. (11) 

i=l 

Considering the string gas as a perfect gas is not only a simplifying assumption but it is also a very 
good approximation. Indeed, the string gas will have diff'erent equations of state depending on the 
temperature. For T ~ Th (the Hagedorn regime), the string gas will behave as pressureless dust 
i.e. P = 0, and the dimensions are nearly frozen at their self-dual radius ~ vc7. The general picture 
is that the universe starts out in a hot dense and isotropic phase close to the Hagedorn temperature, 
then under the effect of some thermal fiuctuations some dimensions will start expanding. The 
number of expanding dimensions is singled by the dynamics. As the temperature drops below the 
Hagedorn temperature (T < Th), the string gas will behave as radiation i.e. P = qE (radiation 
regime). In the following, we will consider two important subcases: the isotropic case where all the 
scale factors are the same, and the anisotropic case where 3 dimensions have already larger radii 
that the remaining 6. 



2.1 The isotropic case 

In this case, the background geometry is described by the metric (^, with ah the scale factors equal 
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ds2 = _dt2 + ^2 g2A(t) ^ ^02^ 0<9,< 27r. 



(12) 



i=l 



2.1.1 Hagedorn regime 



As we said before, at temperatures close to Tjj, the system will behave as pressureless dust. Fur- 
thermore, during this phase, the free energy vanishes so the strings will have a constant energy given 

by 

Eo , (13) 






where S is the entropy. In this regime, and considering the case where all the scale factors are equal, 
one can find an approximate analytic solution for the system of Eqs. (|5l7j) which is given by [5, 14] 
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(^0* + 2B- 2VdA){B + VdA) 



{Eot + 2B + 2VdA){B - VdA) 
A and B are integration constants related to the initial conditions and are given by 



(14) 
(15) 

(16) 



and d is the number of space dimensions (d = 9). This analytic solution is in very good agreement 
with the solutions obtained numerically and is very useful to predict the asymptotic evolution as 
i — > oo. From Eqs. (jTl)) and (fT3|) . we see that a{t) tends to the asymptotic value [14] 
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On the other hand, the shifted dilaton, and thus the dilaton, rolls monotonically to the weak coupling 
regime. 



2.1.2 Radiation regime 

As the temperature decreases below Th, the Universe enters the radiation regime. The free energy 
of the string gas receives two main contributions corresponding to: (i) the string massless modes, 
which we refer to as radiation contribution, and (ii) the string massive modes which we refer to as 
matter contribution. In this paper, we will only consider massless modes since as it was seen in [14] 
they provide the main contribution. The total energy reads [14] 



E[tl = ^D{OfT^'^+^ 
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(18) 



where D{N) is the degeneracy factor (-D(O) = 16) and Cid+ 1) is the Riemannian zeta-function. 



The radiation contribution is characterized by the fact that the d dimensions exhibit a radiation- 
hke equation of state, -P™*^ = E^"''^ /d. By demanding that the universe undergoes an adiabatic 
evolution, 

ls = 0, (19) 

it follows that 

/? = /3o-. (20) 

One expects that at sufficiently late times the dilaton has been already stabilized at weak coupling, 
this means, as noticed first in [5], that Eqs. (|5l7j) admit only solutions with a radiation type equation 
of state. 

Introducing the two variables x = \ and y = ifj m. the case of d = 9, we can write Eqs. ©-© as 

X = xy + -w{y'^ - Qx^) , (21) 

y = l^' + \y'- (22) 

The fixed point of the system corresponds to i = and y = 0. As shown in Ref. [5] the attractor 
solutions satisfy the relation x = —wy. Substituting this relation in Eq. H21|) . we obtain the dif- 
ferential equation x = —(1 + 9'w'^)x'^ /2w. Integrating this equation gives the following asymptotic 
solution: 

In the case of radiation (w = 1/9) one has a oc t^'^. Since the time-derivative of the dilaton 
(j) is given by = {ip + 9A)/2, the dilaton is stabilized for w = 1/9, after the system enters the 
attractor stage characterized by A = —wip. Note that if we consider the d = 3 spatial dimensions 
with w = 1/3 we recover the evolution of the scale factor in a radiation dominated universe [14] 
(a(xti/2). 

Since the condition, 2w < 1 + 9^^, is always satisfied independently of the equation of state, the 
universe does not exhibit an accelerated expansion in the absence of a' corrections. Equation H23|) 
also shows that both A and ip decreases toward zero with time. 

2.2 The anisotropic case (3+6) 

The background geometry in the anisotropic case takes the form 

3 6 

ds^ = -dt^ + -^s X] "?(*) ^^»^ + ^^sYl ^?(*) ^^?' < 0» < 27r and < t?» < 27r, (24) 

i=l 1=1 

where ai{t) = a{t) = e ^ '' is the scale factor of the large dimensions and 6j(t) = b{t) = e^^ '' is the 
scale factor of the small dimensions. Introducing three variables x = X, y = ip and z = /i, we find 
that Eqs. (|5l7j) with an equation of state 



E 



Px = -, P^ = 0, (25) 



where P^ (Px) is the pressure along the small (large) dimensions given by Eq. ®, give 

x = xy- -x^ + -y^ - z^ , (26) 

V = \x^^\v'^'iz\ (27) 

z = yz. (28) 

In Ref. [14] it was shown that the attractor solutions satisfy 

2/ = -3x, z = 0. (29) 

As in the isotropic case, the time-derivative of the dilaton, (/> = (y + 3x + 6z)/2, vanishes along the 
attractor solution. Hence the dilaton is fixed in this regime. 

Substituting Eq. (|^U)) in Eq. (|^H)l , the variable x satisfies the differential equation i, = 2x^ . Then 
in the attractor regime we obtain 
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A = — , /i = 0, a(xt' , 6 = const , tp = —— . (30) 

Hence the late-time evolution for the large dimensions is described by a radiation dominated phase 
with a fixed dilaton. Since the pressure along the small dimensions vanishes, this leads to the 
attractor with /x = const. Hence the large dimensions expand due to the presence of radiation 
whereas the small dimensions can be kept small relative to the large dimensions, provided that the 
initial conditions satisfy Aq > /io and ao > bo. 

If one includes the massive modes they contribute with a pressure for the large and small di- 
mensions. The effect is such that the winding modes lead to a positive (negative) contribution to 
the large (small) dimensions and the momentum modes give a positive contribution to both large 
and small dimensions. However at the self-dual (initial) radius the pressure due to the winding 
modes in the small dimensions compensates the pressure due to the momentum modes resulting in 
a stabilization of the small dimensions. On the other hand the large dimensions keep expanding. 

3 The BV scenario including a' corrections 

Before including higher curvature corrections and studying their effects, let us discuss which parame- 
trization of the a' correction will be used. It is well known that the low-energy dynamics of the 
massless string modes can be described by an effective action. This action is defined as a field-theory 
action generating the S-matrix (scattering amplitudes) which coincides with the massless sector of 
the (tree- level) string S-matrix [38-41]. Since the perturbative S-matrix is invariant under local 
field redefinitions of the background fields there exists a large class of effective actions, which all 
correspond to the same string S-matrix [42,43]. The particular choice of field variables corresponds 
to a particular choice of renormalization-group scheme, for example, "the cj-parametrization" is the 
one most suitable for comparison with the a-model and the "s-parametrization" is more appropriate 
for comparison with the string S-matrix. Such redefinitions affect the higher order in a' terms of 
the target-space effective action. However the equivalence theorem states that there exists a choice 
of "relative" scheme (i.e. of field redefinition) such that the perturbative solutions of the effective 



actions coincide [44-48]. Therefore for simplicity we will choose to work in the cr-parametrization 
scheme. 

We consider the action of type IIB strings in a ((i+l)-dimensional space-time including the 
leading a' corrections, which in the cr-parametrization is 

S = Sq + Sa' + Sm , (31) 

where 8^' is given by [43,49] 



s 1 



'^'^d+1 




(32) 



In type I and heterotic string theory, the leading higher curvature correction is of 0{a') and it is of 
the Gauss-Bonnet type. In type IIB, it turns out that, due to the higher number of supersymmetries, 
the leading higher curvature correction oppears only at 0{a'^). 

Our aim is to find solutions to the equations governing the dynamics including higher curvature 
corrections and to compare them with the corresponding solutions without the corrections. We would 
like to see how diff^erent the solutions are and when the corrections begin to be most relevant. These 
corrections are considered in all cases as perturbations to the original equations and are expected to 
be relevant only at early times where curvature can be large. If this were not the case, one would 
have to include next-to-leading corrections, which is an infinite series in powers of a' . Thus, if next- 
to-leading terms are relevant this means that curvature is so large that the supergravity description 
will break down. 

For our purposes it is important to recall the following two facts. First, since the equations 
including higher curvature corrections are of higher order with respect to the case without the 
corrections and very complicated, it is very difficult to find analytic solutions and thus we will rely 
on numerical integration, which gave excellent results in the standard case. Second, the corrections 
contain higher derivatives and higher powers of derivatives not present in the original equations. This 
fact is important for the structure of the solutions. Let^('^^(t) = {A^ , V }(0 be the solutions of the 
equations without higher order corrections and let '^^'^'{t) = {A^ , il)^'^'^{t) be the solution including 
the corrections. Due to the above, the latter do not reduce to the ones without perturbations when 
taking the limit a' — > 0, i.e. they are not perturbative solutions. This fact makes harder for us to 
find solutions that are physically meaningful. In order to find such solutions, we consider solutions 
that are perturbations of the original ones. We call these solutions perturbative solutions. 

Now, let us outline our strategy for finding the perturbative solutions. Given the above discussion 
about the properties of the equations, we will search for solutions that are perturbations over the 
original ones. Namely we will look for solutions of the form 

^W(t) = ^(0)(t) + a'3 ^(P)(t), (33) 

where ^'^^(t) is the perturbation of the solution ^''^^(f). The consistency of the ansatz Eq. (|33|) 
implies that a'3^^^'(t) <^ ^^'^'(t). Moreover, the ansatz (|33|) guarantees that the solution found 
will reduce to the non-perturbed one as a' — > 0. Plugging this ansatz into the equations and 
expanding all the terms in powers of a', we obtain differential equations that determine ^'■^^(f) 
knowing '^^^>{t). Then, we solve these equations by imposing boundary conditions such that the 
eff'ect of the corrections at late time is negligible i.e. '^^^'{tf) = 0, where tf stands for some late time. 
This corresponds physically to impose that at late times the asymptotic behavior of the solution is 
the same as the unpertubed one: at late times the dynamics must be governed by the R term. 

8 



3.1 The isotropic case 

Let us first study the situation in which the radii of the 9 dimensions are all equal. For the action 
(|31|1 we obtain the following equations of motion: 

- 9A2 + -02 + a'3/i = e'^E, (34) 

A-AV' + a'3/2 = ^e^P, (35) 

V^-9A2 + a'V3 = ^e^S, (36) 
where the correction terms are given by 

/i = C(3) [1890A^ + 216A^V' + 396A5V'A - 594A^A2 + 216A3V'A^ - 144A2 A^ 

+m\i)\^ + 27A^ - 396A^A(^) - 432A3AA(=^) - lOSAA^A^^^^] , (37) 

/2 = 2C(3)[60A^V'-36A^^A + 6A^(V'^-70A-V^) + 22A3A(-3V'A + 4A(3)) 

+4AA2(-5V'A + 12A(^)) + 11A^(V'^A -\ij- 2V'A(3) + a(^)) 

+6A2(V'2a2 + llA^ - A^V' - 4V'AA(3) + 2A(3)2 + 2AA(^)) 

+A(V'^A2 + SA^ - '}?i) - eV'AA^^) + 6A(^)2 + 3XaW)] , (38) 

h = 18C(3)[60A^ + 6A^'(/' + 6A^A-11A^A2-2A2a3 + A^ 

+AA2(V;A - 3A(3)) + Q\^\{'4,\ - 2A(3)) + 11A^(V'A - A(3))] (39) 

In the same way as Eqs. ((SJITI) are invariant under T-duality, Eqs. (|34l36p should be invariant under 
the same symmetry. However the duality transformations Eq. (|llj) are modified as 
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Xi ^ -\i + a'3/(Ai, 0) and </. ^ - ^ A^ + a'^'giX,, 0) , (40) 

where /(Aj,(/)) and g{\i,(f)) are complicated functions which we will not calculate^ . 

Plugging the ansatz (|5n)) into the equations of motion Eqs. (|51l7j) . we obtain differential equations 
for ^(P)(i), that are valid up to 0{a'^). They are given by 

_18A(o)a(p) + 2V'(°)^(p) + /i(AW,V'W) : 
A(P) - (aWV'(p) + a(p)^(°)) + /2(AW, V'W) : 

V;(p)-18AWa(p) + /3(AW,V'^°)) : 

We then integrate numerically these equations, using as an input both the analytic expression and 
the numerical solution for '^^^'{t) when available. 

3.1.1 Hagedorn regime 

As explained in section 2.1.1, during this phase the energy is nearly constant and the pressure is 
zero. We adopted initial conditions around Eq ~ 1000 a'~^' 2, which corresponds to a high energy- 
Hagedorn regime, oq ~ 1 and we have chosen the shifted dilaton to satisfy the condition e^ <C 1 to 

^See e.g. [55,56] for an explicit calculation of the duality transformations for particular metrics. 
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ensure that the string couphng constant is initiahy small and hence the perturbation theory and the 
ideal gas approximation are trustable. We have carried out simulations for a wide variety of initial 
values. Given A^'''(0) and A'-'^^(O) the value of tp^^^O) is fixed by the Hamiltonian constraint ©• 

Figure ^ shows the evolution of the scale factor for the unperturbed solution a^^' (t) and the 
corrected solution a^'^'{t). At late times the effect of the corrections is negligible and the evolution is 
perfectly described by the unperturbed solutions Eqs. (|14I15|) . As one approaches smaller times, the 
perturbed solution starts growing until it dominates completely the unperturbed one. The higher 
curvature terms make the scale factor become larger than in the case without corrections, thus 
causing the expansion rate to become smaller. This is due to the fact that f2{^^^,'P^^') in Eq- l|42|) 
acts as a negative potential at early times, vanishing at late times. This can be understood from 
Eq. H38() : a numerical evaluation shows that at small times the dominant contribution to f2{^^^' , V' ) 
is given by the the third and fourth derivatives of A, resulting in a negative contribution. Under the 
influence of this term, the expansion rate A^'^^ grows, while in the case without the corrections A^'^^ 
decreases, until it reaches the point where the corrections are not relevant anymore. At this point 
A^'^) changes sign and A'-'^^ becomes indistinguishable from A*-*^-*. Therefore the a' corrections act as 
a driving force for the scale factor. On the other hand, for negative initial value of Aq, the effect is 
the opposite; the scale factor is smaller and the expansion rate is bigger, because now /2(A^''\ V ) 
acts as a positive potential (damping force). 

Concerning the dilaton, in the uncorrected case it starts from its initial value at (p^^' = (po, in 
the case of Fig. ^00 — —2.5, and decreases to more negative values. However, when we include the 
corrections the corrected dilaton gets a positive contribution at early times, becoming less negative, 
thus increasing the string coupling, but remaining always in the weak coupling regime. 

It is worth to notice that the effect of the corrections might help to make the WW modes 
annihilate more efficiently. Indeed, including the corrections leads to both a decrease of the expansion 
rate and to an increase of dilaton, and thus the string coupling, which might help to circumvent the 
"all or nothing" problem in [4]. 

How close to t = one can trust the perturbed solution? As one approaches the origin the 
perturbation becomes bigger and therefore our analysis is not valid anymore. Moreover the Hamil- 
tonian constraint must always be satisfied, with the correction terms included as well. It turns out 
that this energy constraint is satisfied up to the inflection point of a^'^'{t). 

Also, the initial conditions play an important role since they determine the magnitude and the 
time of appearance of the effect. The larger the value of A^'^^(O) the larger the correction will be and 
therefore the effect will be noticable at later times. This makes the next-to- leading corrections in 
the a' expansion larger and one would have to take them into account. 



3.1.2 Radiation regime 

In this case we will consider Eqs. (|41II43() with the equation of state for radiation due to a gas of 
strings in a 9 dimensional space, P = gE. In this case the energy is given by Eq. (|18() with d = 9. 

Figure[21 shows the behaviour for the corrected and the uncorrected solutions. As in the Hagedorn 
regime, the late time evolution is perfectly described by the uncorrected equations of motion Eqs. ^ 
I7|) and the effect of the corrected terms is visible only at very small values of time. 
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Figure 1: The evolution of the scale factor a(t) and the expansion rate A(t) for the initial conditions 
A(°)(0) = 10-3, AW(0) = 0.2, V^°HO) = -5 and Eq = 10^ in the Hagedorn regime. 

In contrast with the Hagedorn regime, under the influence of f2{)^^^' ,iIj^^'), the expansion rate 
decreases more rapidly, until the time where the corrections are not relevant anymore. There is 
no inflection point in this case. Once again, this can be understood by inspecting the sign of 
/2(A''^','0 ) in Eq. H42|) which in the radiation regime is positive. Indeed, the expansion rate in 
the radiation regime is bigger than in the Hagedorn regime, and therefore in Eq. ()42() . not only the 
last terms, which involve A'^' and A^^^ will contribute, but also the remaining terms, resulting in a 
positive contribution. For negative values of A''^^(0), the result is the same. 

The dilaton in the uncorrected case starts from an initial negative value, for instance for the 
initial conditions in Figure |2l (/)o — —9.15, then it grows until stabilizing at weak coupling (f) ~ —4. 
The effect of the corrections at early times is to give a negative contribution thus decreasing the 
string coupling constant. 

To summarize, there are two main differences with respect to the Hagedorn case. First, the 
effect of the corrected terms is stronger in the Hagedorn case than in the radiation case, i.e. for 
small values of Aq the effect due to the corrections is noticable in the Hagedorn regime whereas it 
is negligible in the radiation regime. Second, in the Hagedorn regime the effect of the corrections 
results in a driving force for the scale factor, while in the radiation regime, it leads to a damping 
force. Furthermore, the dilaton gets a positive contribution in the Hagedorn regime and a negative 
one during the radiation regime. 



3.1.3 Fixed points 

Even though the exact solutions to Eqs. H34II36() are very difficult to obtain, one can look for the 
presence of fixed points. Introducing new variables A = yi, A = y2) A^^^ = y^ and tp = y^, one can 
rewrite the Eqs. H34II36|) in an autonomous form: 



yi = 2/2 , 
y2 = 2/3 , 



(44) 
(45) 
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Figure 2: The evolution of the scale factor a{t) and the expansion rate A(t) for the initial conditions 
A(°)(0) = 1.33, A(°)(0) = 0.4, V^°nO) = -30 and (3q = 15 in the radiation regime. 



U = \yl + Ivl - a"h + la'^fi 



where 



F 
G 



2C(3y^(llyt + 12yfy2 + 3y2'), 

-y2 + yiy4 + C(3)a '(6y? + Uyfy2 + Gyfyl + yi)(y| + 9y?) 

+f(-9y? + y| + «'Vi)-a'V4, 



and fi is defined in the Appendix. 

The fixed points of this system can be obtained by setting iji = Q {i = 1,2,3,4) [53,54]. 
gives 



9yf + yi = 270C(3)a'-^y^ , 

2/12/4 + |(-92/? + yl) + 6C(3)a'3 [-20y[y4 - 2ylyl + |(315y? + myly^) 



Combining these two equations, we find 

[yi - 3C(3)a'3(40yI + Ayly^)]{yi - 9wyi) = . 

In the rest of this subsection we shall set a'^ = 1 ^. 
When 2/1 - 3C(3)(40y[ + Ayfyi) = 0, one obtains 

(ia) de Sitter : A = 0.8282, ip = -8.1043 , 
(ib) anti de Sitter : A = -0.8282, ip = 8.1043 , 



0. 



(46) 
(47) 

(48) 
(49) 

This 

(50) 

(51) 

(52) 



(53) 
(54) 



3^/3 



a' can be recovered by C(3) -^ a C(3). 
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where we used ^(3) = 1.20206. These fixed points exist for any value of w. The point (ia) is a de 

Sitter solution in which the scale factor grows exponentially (a oc e ), whereas the point (ib) is an 
anti de Sitter solution in which a decreases exponentially. On the other hand, when y^ = 9wyi, we 
obtain the following fixed points: 

l + 9u;2\i/6 . ^ /I + 9^^2x1/6 



(iia) de Sitter : A = ( -^^^ j , ^ = 9w [-^^^ ) , (55) 

l + 9y;2xl/6 /I + 9^2x1/6 



(iib) anti de Sitter : A = — ^, , , ib = —9w ^, , , (56) 

^ ' V 3oc(3) y ' ^ \^ 3oc(3) y ' ^ ' 

which are dependent on w. Solutions with the fixed points (|53II56|) correspond to non-perturbative 
in a' solutions, i.e. {X,ip) ~ a'^^'^. 

On the other hand, let us note that there also exists a perturbative Minkowski fixed point: 

(hi) Minkowski : A = , 1^ = 0. (57) 

In the following we shall study the stability of these points against perturbations and then clarify 
the attractor of the system. 

3.1.4 The stability of fixed points 

( f) 
Consider small perturbations, 6yi, about the fixed points yf derived in the previous subsection, 

i.e. yi = y\ + 5yi. Using the fact that y^ = 2/3 = and G^^' = at the fixed points, we obtain 
the following perturbation equations: 

5yi = 6y2 , (58) 

Sm = Sy3 , (59) 

6y3 = j6G, (60) 

Sm = 9yi6yi + 2/4% - Sfs + 2*^/1 • (61) 

Here 5fi and 6G are given by 

444 
Sfi = ^Ci(5yi , 5/3 = ^ di6yi , 6/4 = "^ eiSyi , (62) 

i=l j=l i=l 

4 

5G = Y,Si5yi, (63) 



i=l 



and Cj, dj, e^ and si are defined in the Appendix. Equation l(HT|l yields 

5y4 = (9yi - lQma?,)a'W)5yi - 108C(3)a'=^y?5y2 + 2/4^^4 • (64) 

Hence the differential equations (|58() - H61|) can be written in the form: 

\yi,y2,y3,m) = M *(yi, 2/2,2/3, 2/4), (65) 
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with 

/ 1 \ 



M 



(66) 



10 

si/F S2/F ss/F S4/F 

\ 9yi - 1080C(3)y[ -108C(3)y6 2/4 J 

If all eigenvalues v of the matrix A4 are negative or have negative real parts, the fixed points are 
regarded to be stable. Meanwhile if at least one of the eigenvalues is positive, the fixed point is not 
a stable attr actor. 

We numerically evaluate the eigenvalues for the equation of state < w < 1 and find the 
following properties for the fixed points (ia), (ib), (iia), (iib) and (iii) introduced in subsection l3.1.31 

• (ia) All eigenvalues are negative. When w = 1/9, for example, one has 

7; = -1.84, -6.26, -8.10, -8.9. (67) 

Hence this de Sitter fixed point is stable. 

• (ib) All eigenvalues are positive. When w = 1/9, for example, one has 

-(; = 1.84,6.26,8.10,8.9. (68) 

Thus this anti de Sitter fixed point is an unstable node. 

• (iia) Only one of the eigenvalues is negative and other eigenvalues are positive or have positive 
real parts. When w = 1/9, for example, one has 

V = 0.28 + l.OOi, 0.28 - l.OOi, -3.34, 3.90. (69) 

Thus this de Sitter fixed point is unstable. 

• (iib) Three of the eigenvalues are negative or have negative real parts, but one of the eigenvalues 
is positive. When w = 1/9, for example, one has 

V = -0.28 + l.OOi, -0.28 - l.OOi, 3.34, -3.90. (70) 

Thus this anti de Sitter fixed point is unstable saddle. 

• (iii) The term F in Eq. (|60() is given hy F = 22(^(3)a'^y| and the term S4 includes the yi term. 
This means that the S4/-F term on the RHS of Eq. H60() exhibits a divergence for yi — > 0. 
Hence the Minkowski fixed point is unstable. 

These results indicate that the de Sitter fixed point (ia) is the only stable attractor. This is 
confirmed by our numerical simulations shown in Figure Ol At this fixed point the time-derivative of 
the dilaton 4> is negative, i.e. (p = {^ + 9A)/2 = —0.3253, which means that the system approaches 
a weak coupling regime. Hence the description of an ideal string gas based on a canonical ensemble 
is valid in this attractor regime. 
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Figure 3: Phase space plot in terms of A and ijj for w = 1/9 with initial conditions A = 0.5 and 
A^^) = 0. The solutions approach the de Sitter fixed point (ia) [A = 0.8282 and ip = —8.1043], which 
is denoted as the point A in the figure. 



We recall that the effect of a' corrections is important when A is order unity (note that we fixed 
a' = 1). The de Sitter fixed point (ia) [A = 0.8282] exists in a high-curvature regime in which the 
a' corrections lead to geometrical inflation. It is interesting to note that the effect of a'^ corrections 
leads to an inflationary solution, while this can not be obtained without a' corrections. This type of 
geometrical inflation often appears in the context of second-order string gravity [54]. Prom the above 
argument the effect of curvature corrections can be important even when the dilaton evolves toward 
the weakly coupled regime. This implies that "the high-curvature regime" is generally different 
from "the strong coupling regime" . However one must be careful because at this value of A other 
subleading order a' corrections may have to be included to perform a complete analysis. 

In the absence of a' corrections the asymptotic behaviour of the scale factor for the radiation 
like equation of state is given by Eq. (PH|) with A and ip decreasing towards 0. Meanwhile we showed 
that the Minkowski fixed point (A, ip) = (0, 0) is unstable in the presence of a' corrections. In fact we 
have numerically confirmed this instability around the fixed point for initial conditions of A, A, X^^\ip 
close to zero. The solutions tend to approach the de Sitter fixed point (ia) provided that A are not 
very small relative to 1. 



3.2 The anisotropic case 

Given the metric ()24() the variation of the action H32|) in the presence of an ideal string gas leads to 
the equations of motion: 



/3i 



3A^ - 6/i^ + V' + a /ii 



e'^E. 



(71) 



15 



\-\iP + a'^h2 = ^e^PA, (72) 

■fi-fi^ + a'^h^ = ^eV-p^, (73) 

V^ - 3A2 - 6/2^ + a'3/i4 = ^e^S, (74) 

where P\ and P^ are the pressure of the string gas along the large and small dimensions, respectively. 
Explicit forms of the hi, /12, /13 and /14 terms are presented in the Appendix. 

As discussed in Sec. 2, the pressures Pa and P^ in the radiation regime are given by Eq. H25() . 
We are interested in the evolution of the large and small dimensions in the presence of a' corrections. 
By plugging in the ansatz (|33|) into Eqs. (|71ll74j) one obtains a set of differential equations analogous 
to Eqs. (|41ll43j) for "^^^>{t). Therefore using the solutions for the uncorrected equations ^''^'(t) we 
can numerically solve for "^^^'{t). 

As it was shown in [14], in the uncorrected case, i.e. hi = 0, i = 1, • • • ,4 the expansion rate 
for the small dimensions /i is exponentially suppressed with the decrease of ifj since in this case 
the pressure P^ is zero and Eq. ()73() integrates to /i = fio^^ ~ ^^- Therefore in the absence of a' 
corrections unless the initial value of |/i| is much larger than unity, the radius b can stay small around 
6 ~ 1. It was also shown that the matter contribution, given by the massive string modes, is such 
that the pressure due to the winding and the momentum modes cancels at values 6 ~ 1 leading to a 
stablization of the small dimensions. Hence we shall only consider initial values of 6 ~ 1. Moreover, 
solutions with growing a and small nearly constant b were found for initial values oq > 60 ~ !> 
tp < and Ao > |/to|- Figures |3 and El show such a solution for initial conditions A('^^(0) = 0.6 and 
il^^'{0) = 0.1. The effect of the higher curvature terms is basically the same as in the isotropic case 
with the same equation of state (radiation). In this case, /12 and /13 in Eqs. (|71ll72jl result in a positive 
potential for A and fi at early times. Thus at initial times the expansion rate of both dimensions 
decreases more rapidly. However, the effect is less dramatic for /i^'^-' than for A'•'^^ This is due to 
the fact that the large dimensions receive also a contribution from the pressure in contrast with the 
small ones where there is no pressure. When A ^ |/i| holds initially it was shown in [14] that in the 
absence of a' corrections it is difficult to keep the small dimensions small relative to the large ones 
and the effect of the corrected terms, at early times, is basically the same as in the previous case. 
For negative values of Aq we get bouncing solutions as in the case without corrections [14]. The large 
dimensions start contracting, then A'*^^ changes sign and then they start to grow. The effect of the 
corrections on the dilaton is basically the same as in the isotropic case with radiation-like equation 
of state. 



3.2.1 Fixed points 

In this subsection we shall derive the fixed points of the system of Eqs. (|71B74p . In analogy to the 
isotropic case the fixed points can be obtained by setting A = 0, A^^^ = 0, A^^^ = and ^ = in 
Eqs. (ini)-(|Zll), which gives 

-3A2-6A2 + V'2^a'3/ii(A = A(3) = A(^)=V^ = 0) = e^^, (75) 

-AV' + a'3/i2(A = A(3) = A(^)=V^ = 0) = ^b^Pa, (76) 
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Figure 4: The evolution of the scale factors a{t) and h{t) for the initial conditions A('^^(0) = 2.6, 
A(o)(0) = 0.6, ^W(O) = 10-3, /i(0)(0) = 0.1, V^°nO) = -20 and /3o = 15 in the radiation regime. 

-AV' + a'=^/i3(A = A(3) = A(4)=V^ = 0) = ^e^P^, (77) 



V'^ + a'^ 



/ii(A = A(3) = A(4)=V^ = 0)-Q'3/i4(A = A(3) = A(4)=V^ = 0) = ^e^S. (78) 



Let us consider the equation of state of the string gas given by Eq. H25() . Then we find from 
Eq. ((77)) that /i = is the solution of Eq. ((77|) . In this case we obtain the following fixed points: 



(A, A, -(A) = (0.84719, 0,0.84719), (-0.84719, 0,-0.84719), 
(1.0852, 0, -4.1566), (-1.0852, 0, 4.1566) . 



(79) 



M\ 



The first and third points correspond to de Sitter solutions for the large dimensions (a ex e ) with 
the small dimensions being stabilized {h = const). If we impose the validity of the weak coupling 
approximation, the third one can be regarded as an ideal solution if it is a stable fixed point. However 
we shall see that this is not a stable attr actor. 

Analogous to the fixed points (|53|) and (|54|) in the isotropic case, there exist two fixed points in 
which both large and small dimensions evolve with the same rate: 



(A, /i, ^) = (0.8282, 0.8282, -8.1043), (-0.8282, -0.8282, 8.1043) 



(80) 



/i in Eqs. (|76|) and (|77|). the RHS of these equations are the same. Since 
3C(3)a'3(40A^ + 4A^^) = for Xif; ^ 0. This is exactly the same equation 



In fact by setting A : 

P^ = 0, we obtain 1 

as the one which we used to derive the fixed points (|53|) and (|54|) . thereby giving the fixed points 

given above. It is worth mentioning that Pa, Py and E vanish in this case. This clearly shows that 

the fixed points ()80|1 exist irrespective of the presence of the string gas. As in section 3.1.3 solutions 

with the above fixed points correspond to non perturbative in a' solutions. On the other hand the 

perturbative Minkowski fixed point (A, /i, ip) = (0, 0, 0) is present as well. 
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Figure 5: The evolution of the expansion rates \{t) and fi{t) for the initial conditions A'''^(0) = 2.6, 
A(o)(0) = 0.6, ^(°)(0) = 10-3, /i(°)(0) = 0.1, V^°H0) = -20 and j3o = 15 in the radiation regime. 

3.2.2 Stability of the fixed points 

From the discussion in Sec. 3.1.4 it is clear that the de Sitter fixed point in Eq. (|8U() is stable whereas 
the anti de Sitter one is unstable. The Minkowski fixed point is not stable as in the isotropic case. 
One can study the stability of fixed points ()79|) by considering linear perturbations. Because of the 
complexity of the a' corrections given in the Appendix, we do not consider perturbation equations 
in the anisotropic case. Instead we shall numerically solve the evolution equations (|71j) - (|74)) for a 
wide range of initial conditions and see the stability of fixed points in the next subsection. It turns 
out that the fixed points (|7n|l are unstable. 

For numerical purpose it is convenient to write the equations in an autonomous form. From 
Eqs. dni) and ((ZH) we get 



V^ = -V'^ + 2^^ + 3A^ - «'^^4 + i^a'^hi . 



Combining Eqs. (f7^ and ((75]) with the use of Eq. (|HT|) . we obtain 



BA(^)+C/i(^)=G, 



and therefore 



A(4) 



2BG - CF 
2^2 - AC ' 



/^ 



(4) 



BF-AG 
2^2 - AC 



(81) 

(82) 
(83) 

(84) 



where 



A 
B 
C 



1 



2;-. 



,•.2 



-C(3)a'3[3(A^ + /i^) + (llA* + 3/2*) + 6(/i^/i + 2A2A) + 3A/i(A/i + 2/i^ + 2/i)] 



■2 



,-.2- 



-C(3)a'^[3A/:t(2(A2 + /^^ + A + /i) + ?,\fj) + 6(/i^A + A^/i + A/i)] 



-C(3)q'3[3(A^ + 5/i^) + 3A* + 55/i* + 6(10/i^/i + A^A) + 2,\ii{\fi + 2A2 + 2A)] , 



(85) 
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Figure 6: Evolution of A, /i and -0 for the initial condition A = 1, A = 0, A^^^ 
0, and ^ = —2. The solution first approaches the fixed point (A,/i, ^) = 



= 0, /I = 0.01, /i = 0, 
11^"^'' = 0, and Tp = —2. ilie solution first approaches ttie tLxed point [A,fi,ip) = (1.0852,0,-4.1566) 
temporaly, but the final attractor is the de Sitter fixed point {X,fi,ip) = (0.8282,0.8282, —8.1043). 



(3) 



and 



G 



-X + \i; + -C(3){A2[6A^ + 3A2/i2 + HA^A + 6X{fi^ + fifl) + 3{2X'^ + /i^ + 2/i^ 
+3fi'^il + /i^A)] + X{X^ + 3/2^ + 3/i^ + 6/u^/i) + 6fi^X{X + /i) + 3A/i^ + ^ilfiXOl + 2A)} 



x{V'^ + 3A2 + 6/i2} + 
1 



-3A^ 



6/i^ + ^"^ + a''^/ii) - a'-^h5 , 



(86) 



-'ll + fiil) + -C(3){/i2[30/i^ + 3AV2 + 55/i^/i + 6/i(A3 + AA) + 3(10/2^ + P + 2A^ 

8 

+3A2A + A^/i)] + /i(5/i2 + 2,'}? + 3A^ + 6A2a) + 6A^/i(A + /i) + 3/lA^ + 3AA/i(A + 2/i)} 
x{V'^ + 3A2 + 6/i2}-a'^/i6. 



(87) 

Here h^ and /le are defined in the Appendix. 

We numerically solve Eqs. (|81j) and (|84|) to understand the attractor behavior of the system. In 
Figure ini we plot one example for the evolution of A, /i and tjj. In this case the initial conditions 
are chosen as {Xi,ili,^i) = (1,0.01, —2). We find that the solution temporarily approaches the third 
fixed point in Eq. ()79p in which only the large dimensions grow exponentially. However the solution 
finally approaches the first fixed point in Eq. ()80j) in which both large and small dimensions exhibit 
accelerated expansion. This suggests that the ideal critical point in which large dimensions grow but 
small dimensions are stabilized is not a stable attractor. In fact we have chosen initial conditions 
close to the third fixed point in Eq. (|79|) and found that solutions repel away from it. This property 
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also applies for other three fixed points in Eq. ((7^1) . 

Before closing this section, let us comment on the consistency of the perturbative solutions that 
we found using the ansatz Eq. (|H5|) and the fixed points analysis above. We found that there exist 
de Sitter fixed points that do not exist when higher curvature terms are neglected. These attractor 
solutions require curvature to be large and correspond to non perturbative solutions in a' . On the 
other hand, at late times our perturbative solutions interpolate to the solutions without corrections 
and ultimately will follow the attractor A = —wip. The reason for that is that our solutions will feel 
the effect of curvature at initial times, then as the universe expands and curvature becomes smaller, 
they relax to the unperturbed solutions. 

4 Conclusions 

We have studied the effect of a' corrections on string gas cosmology. Starting from the low energy ef- 
fective action of type IIB strings including leading order 0{a'^) corrections, in the (T-parametrization, 
we derived the modified equations of motion and studied the adiabatic evolution of the universe. 
With the assumptions of an initially homogeneous toroidal universe and in thermal equilibrium, 
filled with an ideal gas of closed strings at weak coupling, we studied the effect of higher curvature 
corrections finding perturbative (in a'^) solutions and compared them to the uncorrected solutions 
previously studied in e.g. [14]. 

We considered two cases: the isotropic case, where all dimensions are taken to be equal and the 
anisotropic where we consider 3 large and 6 small dimensions. For the isotropic case we considered 
two regimes: a) a high energy (temperature) regime, Hagedorn regime, where the gas of strings is 
described by a pressureless dust equation of state; b) a low temperature regime, radiation regime, 
where the gas of strings is described by a radiation-like equation of state. In the Hagedorn case, 
we found that the correction term acts as a driving force, such that the expansion rate increases 
compared to the uncorrected case. On the other hand, in the radiation regime we found that the 
correction term has the opposite effect, it acts as an initial damping force such that the expansion 
rate decreases faster. We found that in general the effect of the corrections is small and it is only 
noticeable at early times. Moreover it was shown that this effect is stronger in the Hagedorn case 
than in the radiation one. We also showed the existence of fixed points, with and without dependence 
on the equation of state and studied their stability. We found a stable de Sitter fixed point, which 
corresponds to a non perturbative solution, for the radiation-like equation of state that leads to 
geometrical inflation. However one expects that at such large value of A next to leading higher 
curvature terms must be taken into account as well. 

In the anisotropic case, in the radiation regime, we found that the correction terms act as an 
initial damping force for both large and small dimensions initially decreasing their expansion rate 
faster compared to the uncorrected case. Even though the effect is really small, we found that the 
effect is weaker for the small dimensions. This tells us that at least to order 0{a'^) the inclusion 
of higher curvature terms do not alter the stabilization of the small dimensions of the unperturbed 
solution for small values of A. Although the analysis was carried assuming 3 large and 6 small 
dimensions, we expect these features to survive in the more general case of d large and 9 — d 
dimensions. We also analysed the fixed points of the exact system finding one in which the large 
dimensions expand whereas the small dimensions are kept constant but is not a stable attractor. On 
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the other hand, we found a stable de Sitter fixed point which leads to an accelerated expansion of 
both large and small dimensions. Once again, as in the isotropic case, this fixed point corresponds 
to a non perturbative solution. 

Our results show that the inclusion of leading a' corrections affects the dynamics only for rea- 
sonable high curvature initial configurations. In such a situation, the next to leading terms are 
expected to become relevant as well. On the other hand, our perturbative approach is only valid at 
regimes where supergravity description is valid i.e. small curvature. Therefore our work has to be 
understood as a modest qualitative study of the effect of higher curvature in the BV scenario, to 
see in which way the dynamics is affected. There are some open issues that deserve further study. 
For instance, it is interesting to perform a detailed study of the prediction of the dimensionality of 
space-time using our framework [57]. 
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Appendix: Explicit forms of ol corrections in the anisotropic 
case 

In Eqs. (|71I74|) the a' are given by 

hx = ■^C(3) [504A^ + 1260AV^ + 1512 X^ fi^ + 1638AV^ + S02AX^ fi^ + 2772AV^ + 4410/i^ + 144A^V' 
8 L 

+U4:X^fi'^i; + 288X'^fi^ij + 288X^fi^ + lUX^fi^Tp + 720/^^-0 - 216AV^A - lUX^fi^X + 144AV^A 
+216A/1^A + 264^^-0^ + 216X^fi'^tpX + SeoAV^V'A + 2l6Xfi'^ijX - SOGA^A^ - SQGAV^A^ - 216Xfi^X 
+36/P + 144A^V'A^ + 72Xf?ijX^ + 72fi^^'>? - OGA^A^ - 72/^^X3 + 24A^A3 + 18A^ + 216A^/i/i 
+144AV^A - 144A^/i^/i - 2l&>? j/'jl + 216X'^ fi^jil + SeOA^^V'/i + 216A2/:t3'0/i + 1320/i^V'/i - 288X^ fiXil 
-648AV^A/i - 288A/i^A/i + 288AVV'A/i + 288Xf?i)X'il - 216XfiX'^il + 72fL'ip'Pfl + SGA'^/i^ - 216A^/i/i2 
-396AV^/i^ - 1980//"^ /i^ + 72X^tijil'^ + 72X'^fiipil'^ + 720Ai^?/'/i^ - 216A/:iA/i2 + 72X'^X'il^ + lOSA^/i^ 
-72X^il^ - 480/^2/23 + i20/iV'/i^ + 90/i^ - 264A^A(^) - 216A^/i2A(2) - SeoAV^A^^) - 216A/a(2) 
-288A3AA(3) - 144AA2aa(^) - UAfi^XX^^^ - 72XX^X^^'> - 288X^ fi/lX^^^ - 288X^1"^ jlX^^'^ - lUfiXilX^^^ 
-72A/i2A(3) - 216AV/"^^^ - 360AVV^^^ - 216X'^fi^fi^^^ - 1320/iV^^^ - 288A2/iV(3) _ 288^/2^^(3) 
-72/iAV^^^ - 144A3w^^^ - 144A2A/i//(^) - 1440/i^/i//(3) - UAXXjlfi^^^ - 360/i/iV^^^l , (A-1) 

/i2 = ;^C(3) 24A^V' - 33/i'^A - 3A^V'A + '(/'^A^ + 8A^ + Stp'^Xil^ + 9X^fl^ + 15A/i^ + 9^^(5/2^^' - 4V'A 
-10/x/i) + GA^C-^^ - 28A - V^) - A^V^ - SX/l^ - GtpX'^X^^^ - e^/i^A^^) + GAA^^)^ - 12^A/i/x(3) 
+12/iA(^V^^^ + 6A/u(^)^ + 3/i[ - StpXiliX + 3/i) + 8/i2A(^) + 3AV^^^2A/i(2A(^) + 7/i(=^))] 
-3/i3['0A(3A + 10/i) - 2(4/iA(3) + A(2A(3) + 3fi^^^))] + 3A2a(^) + Sjl^X^^^ + 3/t^[V'^A + A(3/i - -0) 
-2'0A(^) + A(^)] + A^[45/i3'0 + IIV'^A - 9/2^ + 3/i2('02 _ 75A - 45/i - V^) - llA^^ - 22V'A(^) 
+3/2(V'/i - 3/i(^)) + llA^'')] + 6A/i/i('') + X^[39fL^ij - 6^(11A2 + /i^) + 6/i^(V'^ - 30A - 26/i - V^) 
+88AA(3) + 6/i/x(3) - 6/i2(2V'(A + /i) + /x^^)) + 6/i(V'^/i - 6A/i - 2/i2 - /iV^ - 2^'/^^^) + /i^''^)] 
+3fi^[X^ + 9A2/i - 4V'/iA(3) + 4A(^V^^^ + 2/iA(^) + A(2V'^/i + 17/2^ - 2/2^^ - 40'/^^^) + 2//(^))] 
+3A2{18/i^?/' + 22A^ + 2A/i2 + 3/i^ + V'^(2A2 + /i^) - 2X'^i; - 'jl^^ + /2^(2V'^ - 39A - 90/i + ^) 
+4A(3)^ - 2/i3(^(3A + 5/i) - /i(^)) + 2/i(3)^ _ 4V'(2AA(3) + /i/x^^)) + /i[ - ^/i(10A + 9/i) + A{2X^J.^^'^ 
+/i(A(^) + 3/|(3)))] + 4AA(^) + 2/i/x(^) + /i2[6/i2 + ^2(A + 3/i) - 3/i'0 - A(6/i + V^) - 2?/'(A(3) + 3/x(3)) 
+A(^) + 3/x('^)]} + A{33/i^V' - 2^^(10X3 + 9A/i2 + 3/i3) + 3/2^(^2 _ 35); _ gg^ - V^) - 3/:(^(4?/'(A + 2/i) 
-3/i(3)) + 3(16A2a(3) + lOA/i/i(3) + /i2(4A(3) + 5^(3))) + f,\- 3-0(3X2 ^ ^g\;^ ^ ^^^2) ^ i2(a(a(3) 
+3//(3)) + /i(3A(3) + 4/x(=^)))] + 6A'[6/i2 + ^'^(x+-il) + A(2/i - V^) - /i^^ - 2^(A(3) + ^,(3)) + A^^) + /x^^)] 
+3/i[6A2/i + 10/i3 + V'2/i(2A + /i) - /i^^; + 4A(3)/x(3) + 2/^(3)^ _ 4^(A/,(3) + /i(A(3) + ^C^))) + 2/iA(^) 
+2/i/x(^) + 2A(9/i2 - /i0 + /iW)]}] , (A-2) 

/i3 = ((3) hlO/i'^V' + 15A^(/iV' - /i) - 180/^^0/1 + 'iip'^X^'il + 15A^/i + 9X2/2^ + 5V'^/i^ + 40/i^ 



4 
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+3X^{9fi'^7p -'ipil + fjiii)"^ - 30A - 18/i - -0)) + 30/i^(V'^ - 49/i - -ip) - S'X^flip - S/i^V' - 12V'A/iA(^) 
+6/iA(3)^ - eV'AV^^^ - 30^/iV^^^ + 12AA(=^V^^^ + SO/i/x^^)^ + fi^{-6ip(X^ + SS/i^) + eAA^^^ 
+440/i//(3)) + /:i(-2^(3A3 + gP/i + 50/i3) + 3(10A/iA(3) + 80/iV^^^ + A2(5A(3) + 4/i(3)))) 
+6A/iA(^) + 3AV^^^ + 15/iV^^^ + 3A'^[13/i^V' + V'^/i + 3A/i + /i^(2'02 - 45A - 39/i - 2^) 
-/iV' + /i(-4V'(2A + /i) + 3A(3)) - 2^fi^^^ + /i(^)] + /^^[-gP + 55(V'^/i - /i-^ - 2V'm(3) ^ ^(4)^^ 
+3A^[30/i^^ - V'/iClOA + 3/i) + 2/2^(^2 _ 26A - 60/i - V^) - 2/i2(^(5A + 3/i) - A^^^) + e/iA^^) 
+8A/i(3) + 4/i/i(3) + 2/i(6p + V'^IA + /i) + A(2/i - i^) - jli^ - 2i;{X^^^ + /i^^)) + A^^) + /x^^))] 
+3X^{33il^ip + 2x/'2A/i + 17A^/i + gA/i^ + /i^ + /i^CV'^ - 90A - 165/i - V') - 2A/i?/; - 2fi^{2tp(X + /i) 
+A(3)) + 4A(^V^^^ - 4V>(/iA(2) + A/u(3)) + A[-V>(11A2 + 18A/i + 3/2^) + 4(/i(3A(3) + ^u^^)) + A(4A(^) 
+3/i(3)))] + 2/iA(^) + 2V(^) + /:i2[6A2 + '02(3A + /i) _ /i^ _ 3A(2/i + ?/;)- 2iP{3X^^^ + /i^^)) 
+3A(^) + ij^% + 3/i2{3A3 + V;2(a2 + io/i2) + A2(2/i - -0) - 4V;(AA(=^) + lO/i/i^^)) + 2AA(^) 
+2(55/2^ - S/i^V^ + A(^)^ + lO/i^^)^ + 10/i;u(^))} - 3A{66/i^A + DV'A^/i + 3V'A/i2 - UA/iA^^^ 
-3/i2A(3) + /i^(-V'A + 3A(3)) - 8AV^^^ - 4A/i;u(3) + /i2[^A(gA + 10/i) - 4(2/iA(3) + A(3A(=^) + ^(3)))] 
-2/23[^2\^ - 2A2 - A(6/i + V^) - 2V'A(^) + A^^^] - ^[lOA^ + tP^X(X + 2/i) + A2(18/i - x/i) 
-4x/-(/iA(3) + A(a(=^) + /x(3))) + 2(A(3)' + 2A(3)^(3) + -fixW) + 2X{3il^ - jmP + A^^) + fi^^^]] , (A-3) 



/i4 = -((3) 24A« + 210/i* + 6A> + 30/i V + A* + 6A''(10/i^ + A) + 30/2^ + 6A^/i^ + 5/i^ 

+/:('^(3A2 - 55/2^) + fi^{-3X^ + 3A2/i - 10/i^) + A^(72/i^ + G/i^-^ + ll^Z-A + 12/i/i - ll^^^)) 
+A^(78/i^ + 12/i3V^ - llA^ - Qfi^(X - 2/i) + 3/2^ + g/i(V^/i - z^^^^)) + 55/i^(V^/i - z^^^^) 
+3A^(48/i^ + 4/i'^V; - 2/:i/i(A + /i) + V'(2A^ + /i^) + A^(V'(3A + 5/i) - 3/i(^) - 5i^(^)) 
-2(2A/i(3) + /izy(3))) + 3/i3(^(A2 + io/i2) _ 2(A/i(3) + lO/ii/^^))) + /i(V'(3A2/i + 5/1^) 
-3(2A/2//(^) + A2i/(3) + 5/i2zy(3))) + A2(132/i^ + Qfi^i^ - 2>? + 6/(2A - /i) + 3A/i2 - 3/2^ 
-3/2^(5X2 + 6A/i + 5/2^) + 3fi^{i){bX + 3/i) - S/x^^) - Sz^^^)) + 3/i(?/'/i(4A + /i) - 2(2Ai/(^) 
+/i(2//(3) + iy(3))))) ^ X{l2fi^X - 6fi^X(X + /i) - 6/iA/i(A + /i) + V'(A^ + 3A/i2) - 9fi\^X 
-/i(3)) - 3(AV^^^ + /^V^^^ + 2A/ii/(3)) + 3/l2(V;A(A + 4/i) - 2(2/i;u(3) + X{fi^^'> + 2i/(3)))))" 

(A-4) 

We checked that these expressions reduce to those in the isotropic case by setting A = //. 
In Eqs. (|86|) and (|87|) the /15 and h^ terms are defined by 

/i5 = C(3);^ 24A^V' - 33/i''A - 3/u^V'A + ip^X^ + 8A^ + 3tp'^Xil^ + DA^/i^ + isA/i^ + 9X^{5fi^ij - AtJjX 

-lOfifl) + 6X^{ip'^ - 28A) - e^Z-A^A^^) - QiPil^X^^^ + GAA^^^^ - 12V'A/i/i(3) + 12/iA(3)^(3) + QX^i^V 
+3/i[ - 3V'A/i(A + 3/i) + 8/i2A(^) + 3AV^^^2A/i(2A(3) + T/i^^^)] - 3/:t^[V'A(3A + 10/i) - 2(4/iA(3) 
+A(2A(3) + 3/i(3)))] + 3/i^[V'^A + 3A/i - 2V'A(3)] + A4[45/i3V' + H^^A - g/i^ + 3/^2(^2 _ ^^-^ 



i{3) 



,(3) 



45/i) - 22V'A^^^ + 3/i(V'/i - 3/i^'^O] + A^[39/i*?/; - 6V'(11A^ + /i^) + 6/i'^(V' - 30A - 26/i) 



+88AA(^) + 6/i/x(3) - 6/i^(2V'(A + /i) + /x^^^ + 6/t(V'^/i - 6A/i - 2/^2 - 2V^^^^^)] + 3/i2[A^ 



/2)) + 6/z(V'^/i - 6A/i - 2/i2 - 2V'^(3))i , o,-,2r\3 
+gP/i - 4V'/iA(3) + 4A(^V^^^ + A(2^^/i + 17/i2 - 4V'^(^^)] + 3A2{18/i^'0 + 22A3 + 2A/i2 + 3/^^ 
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+V'^(2A2 + /i2) + /i4(2V'^ - 39A - 2(45/2)) + 4A(^)^ - 2fi^{i;{3X + 5/i) - ^(^)) + 2^(3)^ - 4V'(2AA(^) 
+/2/i(^)) + /i[ - V'/i(10A + 9/i)4(2V(3) + /2(A(^) + S^u^^)))] + fi^[6il^ + V'^(A + 3/i) - A6/i - 2V'(A(^) 
+3^(^))]} + X{33fi^ip - 2V'(10A^ + QA/i^ + 3/2^) + 3/i^(V'^ - 36A - 66/2) - 3fi^i4ip(X + 2/i) - 3^^^)) 
+3(16A2a(^) + 10A/i/i(^) + /22(4A(^) + 5^^^))) + /t^[ - 3V'(3A2 + 18A/2 + 11/2^) + 12(A(a(3) + 3fi^^^) 
+/2(3A(^) + 4^(3)))] + 6/i3[6/22 + V'2(A + /i) + 2A/2 - 2V'(A(^) + ^(^))] + 3/i[6p/2 + lO/i^ + ip^il{2X 



+/2) + 4A(3)^(3) + 2//(3)' - 4^(Am(=') + /2(A(3) + ;,(3))) + 18A/22]} 



(A-5) 



1 



C(3) 210/i V + 15A''(/i^ - /i) - 180/i^V/i + 3V' A^/2 + 15A-^/i + 9A^/i^ + 5V' /i"^ + 40/2* + 3A^(9/t^V 
^ L 

-V-Zi + ACV"^ - 30A - 18/i)) + 30/f^(V'2 - 49/2) - 12V'A/2A(^) + 6/2A(^)^ - 6V' A^^^^ - 30^-/2^^^^ 

+12AA(^V^^^ + 30/2/x(3)^ + /t3(_6V)(A2 + 55/2^) + 6AA(3) + UOjlfi^^^) + /t(-2V'(3A3 + 9X'^il + 50/2^) 

+3(10A/2A(^) + 80/2V^^^ + A2(5A(^) + 4;u(3)))) + 3X'^[l3il^tp + V'^/i + 3A/2 + A^(2V'^ - 45A - 39/2) 

+/i(-4V'(2A + /2) + 3A(^)) - 2V'/u(^)] + /i'^i-gP + 55(^-^/2 - 2^'^^^))] + 3A^[30/f^V' - V'/i(10A + 3/2) 

+2/f^(V'^ - 26A - 60/i) - 2/i2(V'(5A + 3/2) - A^^^) + 6/2A(^) + SA^^^^ + ijlfi^^^ + 2/t(6A^ + V'^(A + /2) 

+2A/2 - 2V'(A(^) + ;u(^)))] + 3A2{33/i^V' + 2V'^A/2 + 17X^/2 + 9A/22 + jl^ + /t^(V'^ - 90A - 165/2) 

-2/t3(2V>(A + /2) + A(3)) + 4A(3V^^^ - 4V'(/2A(3) + V^^)) + ^[-^-(11X2 + 18A/i + 3/2^) + 4(/2(3A(3) 

+A/(^)) + A(4A(3) + 3/i(3)))] + a2[6A2 + V'^(3A + /i) - 6A/i - 2V'(3A(=') + a^^^^)]} + 3/t2{3A3 + V'^(P 

+10/2^) + 2p/2 - 4V'(AA(^) + 10/2;u(^)) + 2(55/2^ + A^^^^ + lO/i^^^^)} - 3A{66/i^A + 9V'A2/2 + 3V' A/22 

-14A/iA(^) - 3/22a(^) + fi\-ipX + 3A(^)) - 8AV^^^ - 4A/2/i(3) + /t2[^A(9A + 10/2) - 4(2/2A(3) + A(3A(3) 

+H^^^))] - 2/^3 [V;2 A - 2A2 - 6A/i - 2V'A(3)] - il[WX^ + V'^A(A + 2/2) + ISA^/i - 4V'(/2A(^) + A(A(''^) 



+fi^'^)) + 2(A(3)' + 2A(3)^(3)) + 2A3/22]} . 
These are obtained by extracting A^^', fi^'^' and ijj terms from /i2 and /13, respectively. 



(A-6) 



h = 2C(3)[60yIy4-36|/f|/2y4 + 6|/?(y|- 702/2) + 22y3y2(-3|/2y4 + 4^3) 

+%i(y2y4 + II2/2 - 4^22/32/4 + 2y|) + y2(yiy| + 8y| - 62/22/32/4 + 62/i) 
+42/i2/i(-5y2y4 + 122/3) + lly\{y2yl - ^y^VA] 



(A-7) 



ci = 9 C(3)a'^ 1680yS-^^^ + 168^5'''^ ^y^^^ 

C3 = -396C(3)a'3y[-^^5, 

di = 18C(3)a'3 

d3 = -198C(3)a'3yi^^5^ 



480yS-^^^ + 42yS^^6y(/) 



,(/)6„(/) 



,(/)5„(/)2 



ei = 2C(3)a'-^ 4202/}^ ^^y^^ + 3QyY>^yX 
63 



-44C(3)«'3yi^)4y(/)^ 



C2 = 396C(3)a'32/S^)5y(/)^ 
C4 = 216C(3)a'32/S^^^ 



^2 = 18C(3)a'3 

d4 = 108C(3)a'3y(^^^ 
e2 = 2C(3)a'3 

e4 = 2C(3)a'3 



6yS-^^6 + iiy(^)5y(/) 



lly(/)4y(/)2 



36y(/)5y(/) 



602/5^^7 + 12yS-^^6y(/) 



4202/i' 

(A 



(/)6' 
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Si = y4-ei + 36C{3)a'-Y,{yi + 9yf - 270C{3)a'-Yi) + 6C(3)Q'^y?(18yi + 2di - ci) (A-9) 



-9wyi + -^ci 



52 = -1 - 62 + na3)a'-'ytiyi + Ot/f - 270C(3)Q'^y^) + 6C(3)Q'^y?(2d2 - C2) + -C2(A-10) 

53 = _e3 + 6C(3)Q'3y6(2d3-c3) + -C3 (A-11) 

54 = yi-e4 + 6C{^)a'%f{2d4-C4 + 2y4)+wy4 + '^C4 (A-12) 
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